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Abstract 

We treat the two-particle Green's function in the Hubbard model using the recently developed 
r-CPA, a hybrid treatment that applies the coherent-potential approximation (CPA) up to a time r 
related to the inverse of the band width, after which the system is averaged using the virtual-crystal 
approximation (VCA). This model, with suitable approximations, does predict magnetism for a 
modified Stoner criterion. The evaluation of the two-particle propagator in the r-CPA requires 
the solution of the pure CPA, within whose formalism the vertex correction and the weighted 
Green's functions are obtained. The dynamical susceptibility, including the vertex correction and 
the weighted scattering by the residual interaction, is calculated and shows a spin wave spectrum 
in the ferromagnetic regime. 

PACS numbers: 71.10.Fd,75.10.-b,75.40.Gb 
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I. INTRODUCTION 



In the previous paper 1 (to be referred to as I) we developed an improved treatment of 
the magnetic state of the simple Hubbard model, called the r-CPA. This was constructed 
as intermediate between the mean field Stoner model, equivalent to a VCA treatment of the 
energy bands, and the Hubbard's alloy analogy method, which is equivalent to the CPA. 
The single-particle Green's functions, and the properties such as density of states and mag- 
netisation which can be developed from them, are given in I and will be used here without 
further derivations. 

Several important magnetic properties in metallic systems, such as for example the suscep- 
tibility, depend on two-particle propagators. The treatment of these is straightforward in 
VCA type theories but within the CPA method a number of complications arise. In par- 
ticular, because of the correlated scattering of the two particles from a single site, a vertex 
correction has to be included. The theory of such systems in the alloy situation was devel- 
oped by Velicky^ in the 1960's and applied to the calculation of the electrical conduction. 
More recently, Elliott and Schwabe 3 ' 4 have shown how to derive weighted averages of such 
propagators which they applied to the theory of electron-hole pairs combined to form exci- 
tons. 

In the present problem the susceptibility is related to the propagation of a pair of electrons 
of opposite spin, but because electrons of different types see different scattering centres 
these theories have to be generalised to the case of a four-component alloy. In addition, the 
residual interaction between electrons of opposite spin on the same site has to be included 
separately and is the dominant cause of the appearance of spin waves within the model. 
These properties are examined in the pure Hubbard CPA case which, as has been pointed 
out^, has only a paramagnetic phase except in special circumstances. The treatment can be 
carried over into the r-CPA, again requiring a matching of the propagators, for t < r treated 
in the CPA with their vertex corrections, to the VCA solutions appropriate to longer times. 
A full treatment of this situation with further approximations to improve the tractability 
does indeed predict the existence of a broadened spin wave branch in the magnetic response 
functions. 

The equation of motion for the two-particle propagator needed for the susceptibility is de- 
veloped in SecUll Sec JIIII treats the two-particle Green's function in the pure CPA. In Sec lIVI 
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the results are combined in order to construct the solution in the r-CPA. The dynamical 
susceptibility is then calculated in SecfVl 



II. TWO-PARTICLE GREEN'S FUNCTIONS IN THE HUBBARD MODEL 

As was done in I we consider the Hubbard model 

H = H + UY^ 4^4i c ih H o = Yl €kC k* c k* (1) 

i k 

for N spins at temperature T = 0. A particular property we want to study in the Hubbard 
model is the response of the system to the application of an external magnetic field. In 
the linear response theory^ the induced moment is given in term of the susceptibility tensor 
X Mi/ (x — x', £ — t'), /i, z/ G {x, y,z}. The susceptibility tensor is completely determined by 
the longitudinal susceptibility \ zz an d the transverse susceptibilities \~ + an d X + ~ ■ We will 
concentrate in this paper on the transverse susceptibilities, which are given as^ 

x -+(x - x', t - If) = i2*0(t - 0<[*~(x, t), a + (x', if)]) (2) 
X +-( X _ x ', t - t >) = i27TB(t - t')([^ + (x, t), If)]) (3) 

where a + (x) = cic^ and cr~(x) = 4i°x\- 

The general Green's function that we want to evaluate are therefore of the type 

Ki jk i(t) ^ {(c+^c^c+c^)) (4) 

The transverse susceptibilities will be found as x ^~( x )t) — — Kxxoo{t) with a similar expres- 
sion for x + ~{ x i t) given by exchanging the f and | in 
The equation of motion for K^i in the Hubbard Hamiltonian reads 

EKijki = Sjk(cti c h) ~ dn( c ki c jl) 

"I - ^ tjn^inkl ^ ^niK n jkl 
n n 

+ U ((( c ti c 3^ji - ^T c J c it) ; 4i c h))e (5) 

where n ia is the usual particle number operator 4 a Ci a . A usual procedure to close the 
equation of motion is to substitute the operators n ia with numbers. In the simplest approxi- 
mation the average number of spins N a is taken. This leads to the VCA treatment (or Stoner 



model). If on the other hand the alloy analogy is applied and the operators are replaced by 
a random number, or 1, we will get the CPA result for the two-particle Green's function. 
However, replacing the operators by numbers must be done carefully. We note that the last 
term in (j5J) disappears when i = j, as both c^Ciic^Cii and c^c^c^c^ independently go to 
zero. This feature is lost if only the substitution n ia — > rii a is made, where n ia is either the 
average or or 1 at random. The equation of motion therefore goes over to 

EK ijkl = 5 ife (cJc Zi ) - ^(4t c jt) 

n n 

+ U (((n jl - mi) cJc iT ; 4 T c u )) E 

- US ij (((n jl - mi) c+c iT ; 4 T c u )) E (6) 

The term artificially introduced is thus simply subtracted again when i = j. 
The evaluation of the two-particle propagator in the Hubbard model using the r-CPA meth- 
ods require the solution of the pure CPA and VCA case. We will first treat the CPA case. 
The VCA case can be easily determined and will be discussed later. 



III. TWO-PARTICLE GREEN'S FUNCTION IN THE PURE CPA 

In this section the equation of motion (jUJ) is treated in the alloy analogy. The numbers 
in © are now either 1 or 0, depending on whether there is a spin a at i or not. The 
notation G a will be used throughout this section instead of G aCPA . 



A. The four- component alloy 

The alloy analogy in the two-particle case is based on the following observation: assuming 
that the last term in © can be neglected, the same equation of motion is obtained if the 
following reduced, quadratic Hamiltonian is used instead of the Hubbard Hamiltonian: 

Hred = l ii c t c io + Yl ei ° c t c i° ( 7 ) 

i,j,a i,<T 

with €ji = rijJJ and = n^U (8) 
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The residual interaction, the last term of will be included in the calculations at a later 
stage. The rest of this section is dedicated to finding an approximation for K^ kl , where K^ kl 
is a solution of the reduced equation of motion 

EK ijki = M C J C U> - 5 u{ c t} c i\) 

+ ^ ] tjnK ink i — ^ ] t n iK n j kl 
n n 

+ U (n jl - rii T ) ((cJc iT ; c+ c u )) B (9) 

The fact that a quadratic Hamiltonian (J7J) is now used presents a number of advantages. For 
instance, it will be possible to use Wick's theorem to decompose the two-particle Green's 
function into pairs of one-particle Green's functions. Moreover, it is now possible to apply 
the CPA formalism. The type of random energies defined in the reduced Hamiltonian (JZJ) 
suggest the following extension of the binary alloy analogy to the two-particle case. As seen 
from (jHJ), an up-spin travelling on the lattice is interacting only on the sites where there is a 
down-spin. This can be interpreted in the CPA as a particle, for example a spin up, moving 
in the lattice and meeting either a "host" (site with no spin down) or an "impurity" (site 
with a spin down). Since a travelling spin up does not see the same 'hosts' and 'impurities' 
as a spin down, we need to define the lattice as being made of four types of sites rather 
than just two. Let n + ,n~,n ± and be respectively, for each site, the probability of finding 
a single electron of spin up, a single electron of spin down, both one spin up and one spin 
down, and no electron at all. The set {n + , n~ , n ± , n^} defines the four-component alloy for 
the system. Thus = n + + is the concentration of "impurities" seen by a down-spin 
and = n~ + n the one seen by an up-spin. The CPA equation can be written in the 
form& 

N-°(U-E°) _ (1 - N-°)E° = 
1 - (U - S ff )F CT 1 + Y!i a F a 1 ' 

where F a represents the trace of G° ' . Each equation (one for f and one for j) is to be solved 

self-consistently for F a and S°". The two equations |T0|) are not independent, as iV^ and 

must be consistent together with the chemical potential for the system to be in equilibrium. 

The quadratic form of (J7J) allows K® ikl to be decomposed into one-particle Green's functions 
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and correlation functions: 



(( c J(*) c it(*); c feT c u)) = 

«(*)Qi)((ciT(*);4 T )) + (4^(t))((4i(t)^n)) (ii) 

Using the relation between time correlation functions and Green's functions, (c^(£)qj_) and 
(ci.Cj|(t)) can be expressed in terms of and G\ at T = 

(c+ CT (0O = ^ due*"* f-~I«i«<w; c+ CT »-) (12) 
i^-jfcj can now be written in terms of the one-particle Green's functions 
K? jkl (E) = — dutf{j,k',u + E)TmG%v,u) 

+— r duIiiG t (j , l jfe;w)G* l (l l i;w-£) (13) 

./-oo 

It is convenient to express the imaginary part of any retarded Green's function G R (E) as 
the difference over the branch cut: 

lim (G {E + ie) — G(E — it)) = 2 i MG R {E) (14) 

(fT3"|) becomes 

K% M {E) = ^f_ M&U, k;uj + E) G% i; w) 

-G%£;;cu + £) CT^^u;) 
+ G\],k;u) G* l (l,i;uj - E) 

-G*1(j,k;u)G*(l,i;u-E)] (15) 

This expression for Kf- kl must now be averaged over all possible configurations of spins. This 
involves the evaluation of quantities such as (G 1 \zl)G^ \z2)) . As the particles are subject to 
statistical correlations on the lattice, (G^ (zl)G^(z2)) ^ (G^ (zl))(G^(z2)). The correction 
to (G^(zl)){G^(z2)) is expected to be small, as the up-spins do not propagate on the same 
lattice of impurities as the down-spins. 



B. Vertex corrections for the four- component alloy 

In this part the vertex correction is developed for a four-component alloy that is con- 
sistent with the CPA. This is a straightforward generalisation of the binary alloy case ob- 
tained by Velicky^.The vertex corrections arise from the statistical correlations felt by the 
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two spins as they move on the same lattice. We adopt Velicky's procedure and look into 
evaluating G^(z 1 ,z 2 ) :— (G 1, (zi)CG^(z 2 )}, where C is a general operator. It will be later 
restricted to diagonal operators, since for the magnetic susceptibility only quantities of the 
type (G\ m {z 1 )G\ nl {z 2 )) are needed. The result, after directly following Velicky's algebra, 
can be summarised as follows 

{Zl,Z 2 )-- , y , — v UDJ 

1 - A{zi, z 2 )A k {z 1 , z 2 ) 

where 

a k ( Zl ,z 2 ) = J2e~ ikR ™ (m\G\z 1 )CG i (z 2 ) \m) (17) 

m 

A k (z l7 z 2 ) = e- tkR -Fi( Zl )Fl(z 2 ) (18) 

m 

F n _ m :=(n\G\m) (19) 
can be written in products of the single-particle functions. The vertex correction 

AM = MWW ,20) 

gives the correlated scattering of the two particles for the same sites. We recall that the 
CPA result for the t n 's is t n (z) = (e„ — E(z)) [1 — [e n — H(z)]F(z)]~ 1 . If C is now restricted 
to a diagonal operator, A k becomes equivalent to a k , and they represent the free-particle 
solution 

= jfYiGUMGfa) (2i) 
q 

The evaluation of A for the four- component alloy follows again on the same lines as for 
the binary alloy in the single-site representation. The average (tl l (zi)t\ l (z 2 )) is obtained by 
summing the contributions from the four species on the lattice. To simplify the notation a 
little the energies z\ and z 2 have been omitted again: 



n 



(-ET) (U-^) 



+ n 



± 



[1 + £TFT][1- (U-El)Fl] 



+ n 



(C/-ET) (-E1) 



[1 - (C/-ET)FT][l + EiFi] 

+ „o (-*)(-*) 



[l-(C/-ET)FT][l-([/-Ei)Fi] [1 + E?FT] [1 + E^] 



(22) 
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This expression can be simplified by substituting the CPA result F a = 
(E CT - N-°U) (U - S' 7 )]- 1 in Defining 

n ■= n ± - N + N~ ri := — (23) 

' N+N-{1- N+){1- N-) 1 ; 

we find 

, t lx E^ (U-^) (U-&) , 

(tlA) = * ^ L V (24) 

Putting ()24j) into (f2"U|) leads to the vertex correction 

a _ W(U-tf)(U-&)n 

U 2 WNin° + fi(Em - N+U^ - N~UEl) 1 ' 

We note for later reference that the vertex correction can be put in the form 

!_ 1 fN+N-n n + N-(l-N+) 

~ n± - N+N- V StS^ + ST([/- Si) 
n-iV+(l-iV-) ^(l-iV+)(l-iV-) \ 
Ei(C/-ET) + (C/-Ei)(C/-ET) ^ 1 J 

In the limit U -»• 0, (44) = (*T) (fi) = 0, so that A tends to zero, n 1 * 1 takes the same value 

as in the purely random case where the two particles up and down are not interacting with 

anything, that is — > N + N~ . In this respect we note that the vertex correction, with its 

factor (n^ — N + N~) in the numerator, gives a measure of the deviation from the random 

value. 

In the limit U — ► oo, E CT — > (-A^< T )[F CT ]- 1 , so that A -> n[n° F^' 1 . In this limit 
the interaction should prevent any double occupation of the same site, so that n — > 0. 
The resulting A is proportional to both the concentration of species f and j, and inversely 
proportional to the proportion of vacant sites. 



1. Effect in the large U limit 

The large U limit is conveniently obtained in conjunction with the limit w — > 0, where 
w is half the band width. This limit is to be treated with caution when dealing with the 
susceptibility, since no transitions are possible if the bands have no width. It is nonetheless 
possible to get an idea of what will be the respective weights associated with the bands in 
the split-band limit. In this limit, the behaviour of the Green's function in the CPA is given 
by£ 

1 - N~ a N~ a 
G ^) = _^ + _ (27) 
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The two-particle Green's function Gq (E\,E 2 ) = Gfi(E\)G^(E 2 ) becomes, in those limits 



r®(F Nmi ■ N] ( l ~ Nl ) | q-N^m (l-iVT)(l-iVi) ^ 

G ° 1 U 2) ~ (Eh - U)(E 2 -U) + E X {E 2 -U) + E 2 (E ± - U) + E X E 2 W 

The transverse susceptibility x" + consists of transitions from the down-spin states to the up- 
spin states. Though (|2*5|) suggests that four peaks will appear, the features of x~ + strongly 
depend on the filling of the bands. In a less than half-filled band, the chemical potential 
will typically fall within the main sub-band. After the energy convolution on (|28jl. only the 
two last terms in (|2"%|) will remain. 

The inclusion of the vertex correction (|2*5j) has the effect of modifying the weights in the 
following manner: 



1 u 2} (E x -U){E 2 -Uy E X {E 2 -Uy E 2 {Ei-Uy E X E 2 1 ; 

Since is always smaller than or equal to N*N^, the first and last terms of (|29*jl are reduced 
compared to (|28|). On the contrary, the two middle terms are enhanced. It is therefore to 
be expected that the vertex corrections will favour transitions from the main band to the 
impurity band and reduce transitions between the impurity bands (if they were allowed) or 
between the main bands. We note that the four weights in (|2*9*j) are actually the respective 
concentrations of species. 



C. Weighted Green's functions 

The concept of weighted Green's functions becomes necessary when different types of 
objects, such as the A and B atoms for the binary alloy, are on the lattice. Weighting in 
the CPA means that we choose from which type of site the particle is starting the motion, 
as well as in some cases the type of site where the motion ends. This results in a statistical 
weight given for each component. The properly weighted Green's function then allows us 
to identify and analyse the effects of each component of the alloy. Recently, Elliott and 
Schwabe^ calculated all the necessary weights for the two-particle propagator in the binary 
alloy. In the subsequent paper—, they showed how to use the weights to obtain a more refined 
treatment of the interacting exciton in this context. The two-particle weighted Green's 
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functions in the four component alloy will be required when the residual interaction in © 
is included in the calculations. Moreover, it will be shown that the one-particle weighted 
Green's functions can be used to determine n ± . 



1. Weighted one-particle Green's functions 

We consider first the one-particle propagator and assume that a particle propagates be- 
tween "impurity" sites i and "host" sites h (as it does in the four- component alloy as well). 
The singly-weighted Green's functions are defined such that 

G(L m) if 1 is an impurity 
G io (l,m) = { (30) 

otherwise 

We also define G m (l,m) as corresponding to the original propagator G(l,m) only if m is an 
impurity, and zero otherwise, and the doubly- weighted Green's function G n (l, m) that is not 
zero only if both I and m are impurities. We define similarly G ho (l, m), G oh (l, m), G lh (l, m), 
G hl (l, m) and G hh (l, m). The weighted Green's function in an alloy of the type H = H + V 
like can be obtained simply by applying the diagonal interaction operator: 

G io = jjG, G oi = Gjj, G li = jpVGV (31) 

An estimation of the averages such as G to is needed in order to study how, on average, 
one-particle properties are made of processes on the different species of the alloy. In the 
CPA the one-particle single weights are found to be^ 

Qio = Qoi = _Q Qho = Qoh = (]_ _ }Q (33) 

G to offers a method for obtaining the average number of sites that are doubly occupied in 
the Hubbard model for the paramagnetic case. At zero temperature, the average number 
of sites per atom with an j spin is given by = (cicjf) = J^ oo dE(—l/7r)'MGl i (E). The 
Green's function G*(E) is evaluated in the CPA. This formulation counts all sites where 
an I spin is present, including sites where both up and down spins are found. Recalling 
the alloy analogy, the f spin sees the sites where a j spin is present as defect sites i, where 
the particle interacts with energy U . The weighted Green's function G\l° is thus on sites 
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FIG. 1: n ± for various values of the particle concentration c in the CPA treatment of the Hubbard 
model. The state of the system is paramagnetic. 

where there is no j spin, filtering out the singly occupied sites. Hence the number of sites 
per atoms n where two spins are found, is, on average and in the CPA approximation 

oo 



n ± = J dE l~ \ TmGi;°(E) (33) 
= y dE (~j lm[F\E)- 



(U -tf(E))Fl(Ey 

In the limit U — > 0, — > N' N^, as expected. The results for in the CPA treatment and 
for the paramagnetic case, = A^, are given in FigO For less than half- filling, that is 
for c = = A^ < 0.5, n decreases rapidly from the random value c 2 to 0. At U = 1, n 
becomes almost negligible. Beyond half-filling, saturation forces a minimum number of up- 
and down-spins to cohabit even at large U, reaching a minimum value of 2c — 1 at U — > oo. 



2. Singly weighted two-particle Green's functions 

Weighted two-particle Green's functions are defined identically to the one-particle Green's 
function. The aim is to estimate (G^" 1 ™ 2 ^™ 2 ™ 1 ), where n\ and n 2 restrict or not the type 
of site % and j can be. The result will be naturally linked to the vertex correction. Elliott 
and Schwabe^ presented the calculations to obtain singly weighted and doubly weighted 
two-particle Green's functions for a binary alloy with clearly defined and fixed impurities i 
and hosts h. They considered the possible weighted Green's functions generated by (GCG), 
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where C was specialised to diagonal operators or sums of diagonal operators. For the 
diagonal operator they obtained (G to CG m ) and (G u CG n ), from which they could deduce 
(G ho CG oh ) and (G hh CG hh ) . Only the following weighted Green's functions will be required 
for the susceptibility: (G]fG$ + ), (G$~G$~), (Gjf n+ G\f n+ ), (G]f n ~ G^~ n ~). For 
completeness, the other singly and doubly weighted Green's functions will also be calculated. 



The calculation of (G Tn °CG lon ) corresponds to calculating (G^ io CG l °' 1 }, as rv^ is acting; 
as an impurity for both the up-spin and the down-spin. The calculation is essentially the 
same to the one made by Elliott and Schwabe 3 and will thus not be repeated here. The 
result is 

(G^ n °CG lon ) = (G^CG 1 ) J 2 (34) 

with D := 1 + S T F T + The calculation for (Gfl n ~°CG ian ~) is carried through in detail 

in appendix [X] as it differs from the Schwabe cases. The other weights are obtained in the 
same way. The results are 

(G 1no CG lon ) = (G Un CG ino ) = (G r CG l ) C (35) 

where n G {n ± ,n + ,n~ ,n } and the following weights were calculated: 



(U - ST) (U-Ei) n± - A/W 

A n-N\l-N l ) 
(£/-£T)£l n±-AW 

A n+(l-A^ T )ATl 
£T([/_£l) n ± — ATT A/4 



r ■= — n0N ' Ni (36) 



5. Doubly weighted two-particle Green's functions 

Elliott and Schwabe 3 obtained the following result for the doubly weighted propagator in 
the binary alloy, with C a diagonal operator: 

/vTvl i a n\ 2 

(G^CG lu ) = (G^CG 1 ) ( J 2 ) (37) 

That is, the weight of the doubly-weighted Green's function is the single weight squared. 
This result is not valid for general non-diagonal operators C, as Elliott and Schwabe showed. 
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The result (|37p also corresponds to the four-component alloy case (G Tn n CG in n ). The 
case (G^ n n CG^ n n ) is treated in detail in appendix [BJ The result can be generalised to 
the other doubly weighted Green's functions. We get: 

(G Tnin2 CG in2ni ) = (G^CG 1 ) i ni i n2 (38) 

where rii,n 2 € {n^, n + , n~, n }. 



4- Remarks about the weights 

The relation between the weights and the vertex correction becomes transparent when 
looking at (J26j) . The four terms summed in A -1 appear to be closely related to the weights, 
with flUJ) rewritten as A~ x = A" 1 ^ + £ n+ + C + £ n± ) and the weights summing up to 
1. 

In the small U limit, £ CT — > N~ a U. The vertex correction (|25|) thus simplifies to U 2 (n — 
WN^). Introducing these approximations into the weights lead to 

C l± -> n ± C~ -> n~ C + -> n + ^ -> n (39) 

In the large U limit the weights retain their energy dependence and estimations are of little 
practical use. It is however interesting to note that using the large U limits S " — > A^ -0 "^ "] -1 
and A — > (n ± - N + N~) [n ^^] -1 , we get (F T and are bounded) 

,n± ^(1-^(1-^) . 

? n (F^U + N^)(F^U + N^) ? 

^0 (40) 
? n (F-£/ + A^) 1 J 

where a G {+,— } or {f, J.} accordingly. In the infinitely large U limit the weight is thus 
solely carried by the contribution from the empty sites. Contributions from occupied sites 
are eliminated as the repulsive interaction grows, preventing encounters of the particles on 
those sites when N P < 1. 



D. The residual interaction 



In this subsection we turn our attention to the last term in (JBJ). So far we have found 
the solution K^ kl of the reduced equation of motion. The simplest way of treating the 
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term —U5ij(({riji — n^) ctcjf, ctcjj,)) e is to average the rija's to N a, s, so that it becomes 
— UmSij((ctcjf, c^Cii))e, with m = N^- — . This is the VCA treatment of this term. 
Calling Kijki the solution of the whole equation of motion (JUJ), it can be checked that K xxx i x i 
can be expressed in terms of the uncorrected solution by 

K xxx > x >(t) = K° xxx , x ,(t) 

/oo 
dt'Kl xm {t-t')K nfKC , x ,{t') (41) 
-°° 

Fourier transforming both in time and space, we get 

-irKiy 

This result, however, is not consistent with the one-particle Green's functions obtained with 
the CPA, as fl4*2j) predicts an instability in the ground state at U = —1/K°(0,0). Since the 
CPA does not allow for ferromagnetic solutions, a treatment of the extra term is required 
that is consistent with the original method. 

A more refined treatment can be obtained when ^ if the extra term is treated using 
the weighted Green's functions developed in the CPA. Instead of taking their averages, the 
operators rii a are again replaced by or 1 at random. The integral equation including the 
extra term now reads 

K xx ,(t) = K° xx ,(t) 

- - E ^ - n ^ / dt ' K ^ - 0^(0 (43) 

where to simplify the notation, we have written K xx , instead of K xxx , x ,. Transformed in the 
energy space, it can be developed into 

K XX ,{E) = K° XX ,(E) K °^ E ) ( n n - "«T) K nA E ) 

fit _ 

n 

+ (^) _E K ^(E) (n n - 4 - n n - lT ) K- in - 2 (E) 

x(n n - 2l -n n - 2l )K° 2X ,(E) + ... (44) 

The terms (n„j — n^) can take three possible values, 0, 1 or — 1, depending on whether the 
site n is an empty or doubly occupied site, a single { site, or a single j site, respectively. 
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The summations over n, ni, n 2 etc in f)44j) can thus be restricted to singly-occupied sites. 
For example, the first summation in (j44j) becomes 



n 

= ~11 K ^E)KIAE) + ^Y. K *n(E)KlAE) (45) 

The restricted summation of full Green's functions can be replaced by a summation over all 
sites of the weighted Green's functions. For example, 

--J2 k °^e)kUe) 

m z — ' 

nej, 

= -^E K ^ n ~(E)K^(E) (46) 

n 

where the summation has to be renormalised by the impurity concentration, as the aver- 
aged weighted Green's function is the total counted on all impurity sites. We recall that 
K®® n (E) represents an energy convolution over G® 7 ^ G% x °- Calling v\ = U/{mn + ) and 
f| = —U/(mn~), the expansion (j4^j) becomes, after Fourier transforming 

K{q, E) = K°(q, E) + J] v s K°> os (q, E)K°> so (q, E) 

se{U} 

+ V s^s- 2 K Q '° S1 (q, E)K°' SlS2 (q, E)K°^°(q, E) 

si,«2e{T,i} 

+ ... (47) 

Such an expansion is very similar to what is done in the case of the excitonic absorption 
in a binary alloy 4 . In this paper, Elliott and Schwabe calculate the optical absorption of 
an exciton, taking into account the interaction between the hole and the electron. This is 
done by writing a weighted scattering expansion comparable to ([47)1 . Similarly to RefJJ, 
()47|) can be expressed in a compact matrix form. The following matrices are defined, where 
all weighted Green's functions are understood to depend on q and E: 



K° 



R 0n+\ 






1 1 


y K 0n- j 






f K 0n+ 


K° 


= 






> K 0n+ 






u 



(48) 
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Since the disorder average has been taken, the following simplification can be made: 

K 0n+n- = K 0n-n+^ and R 0n = R 0on _ R 0no for any n gcattering expans i on gTj) 

becomes 

K(q, E) = K°(q, E) + K^UK + K^UKWK + ... 

CO 



= K°(q,E) + K° uY^[K°U) n K° (49) 

n=0 

As f!49|) represents a geometric series, the equation takes the elegant form 

K(q, E) = K°(q, E) + U[l - i^°f/] _1 if (50) 



The result can be written in terms of the elements of the matrices defined in ()48|) . so that 
(|5T7|) goes over to 



(51) 



K(q,E) = K°(q,E)+ 

vi(K 0n ~) 2 + t» T (K°" + ) 2 + vfViKi 
1 - viK Qn ~ n ~ - v^K 0n+n+ + v^iK 2 

with the following definitions: 

j£ = 2K 0n+ K 0n ~ x 0n+n ~ 

— ^x 0n ') 2 K 0n+n+ — (^x an+ ) 2 K ann ~ (52) 
K 2 = K 0n+n+ K 0n ~ n ~ - [K 0n+n ~f (53) 

This expression can be greatly simplified if the small U limit is taken. In that case, the 
f/ 2 -terms are neglected, and the weights are energy-independent (cf Using the property 
(jBSj) for the double-weights, (jKTjl becomes 

„ ^ (54, 



so that the standard result (|42|) is retrieved. 

In the limit U — > oo, we recall that the weights for n + , n - , as well as for n ± are all zero. 
However, the weights in this case are not energy-independent and cannot be taken out of 
the convolution represented by K 0n+n+ (q, E) or K 0n n (q, E), and the behaviour is difficult 
to predict. 
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E. The paramagnetic case 



As the CPA does not allow for ferromagnetic ground states, the paramagnetic case where 
n + = n~, = and m = must be examined separately. In this case, the aver- 
aged weighted Green's functions with respect to n + and n~ are expected to be identical. 
However, the condition that K 0on+ is equivalent to K 0on is not fulfilled by the expres- 
sions for the weights as they are given by ([36)1 . This is due to the fact that the denom- 
inator of G 2 ' n ~ oon ~ (zi, z 2 ) contains (U — S^(^ 2 ), whereas G 2 ' n+oon+ (zi, z 2 ) contains 
T$(zi) (U — (z 2 )) ■ The energy convolution therefore produces two different results for 
K 0n+ and K 0n , in spite of the physical symmetry in the problem. The symmetry is how- 
ever restored if K°_ + is considered together with K9_. In the paramagnetic regime we have 
^-0n+ _ /fQn , with similar equivalences for the doubly-weighted Green's functions. It is 
thus possible in principle to obtain the limit of K xx (q, E) := K_ + (q, E) + K + _(q, E) when 
m — > 0. To simplify the notation, we abandon the wave vector and energy dependences. 
The first term of the expansion of K xx from (|47J) . K\, is given by 

Ki= u( {Kl + + f _ (KT + ) 2 + (Kf_) 2 _ (Kl_)A 
m I n + n~ n + n~ J 

The two terms in the middle and the two remaining terms on the right-hand side of (|55j) 
become identical in the paramagnetic limit and can be paired off. We need only consider 
one pair, for example lim m ^o = U/m (^K^) 2 / n+ — (K\S) 2 1 n J ■ We define iV so that 
iV^ = N — m/2 and N^- = N + m/2. Similarly, n + = n — m/2 and n~ = h + m/2. In the 
paramag netic limit K_, = K"_ -> K n . Both K_, and can be developed in powers of 



m so that 

(J 

n + h dn \ h J dm 

The two terms together give 



n + n \ n 2 n dn 

Due to the self-consistent nature of the CPA, the quantity dK n /dh is very difficult to 
obtain. The self-energy varies in a complex way with the concentration of impurities, as 
the background itself is affected by a small modification in the number of impurities. The 
weights themselves are also dependent on the self-energy and the concentration of species 
in a complex manner. In two cases, however, the quantity dK n /dn can be established. In 
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the limit U — > it is easy to see that dK n /dn — > K n jh. The same result is obtained in the 
split-band limit, since in this case K n (z\, Zq) is given simply by nj\z\{zi — U)]. In the limit 
of small and large U, the first term becomes 



It is easy to see that the subsequent terms in U can all be paired in the same way. Extrap- 
olating the result (|58|) to the higher order terms, we find 



Numerical results are given in the next subsection. 
F. Numerical illustrations 

The numerical applications of the formalism developed for the CPA susceptibility are 
presented here. We use a three-dimensional, simple cubic band at T = 0. We first evaluate 
the susceptibility Xo x (q,E) without taking into account the residual interaction in order to 
illustrate the effects of the vertex correction. The results x xx (q> E) that include the residual 
interaction are discussed next. 

The susceptibility is XQ X (q,E) = Xo + (<L E) + Xq~{q.iE) which, in the paramagnetic 
ground state of the CPA, becomes simply x xx (<l,E) = 2xo + (q, E). Xo + ( c l>E) is obtained 
as Xo + {<liE) = —K°(q,E), where K°(q,E) is calculated from the combination of Green's 
functions l)15J). Each of these Green's functions is obtained from the CPA solution (fTH|l 
with the vertex correction (|25j) . The quantity n is determined from the weighted Green's 
functions (}3*3*j) . The transitions involved in calculating the Green's function K°(q,E), at 
q = and without the inclusion of the vertex corrections, are illustrated in Fig|21 It shows 
the density of states in a split-band case, U = 3.0, iV + = N~ = 0.35 and /i = 0.02. In 
this paramagnetic case the up-spin band is identical to the down-spin band. Transitions 
occur between the occupied states (shaded area in FigEl) of the down-spins and the empty 
states of the up-spins. As the chemical potential falls within the lower sub-band A, 
transitions start at zero energy E. A first peak in the imaginary part of Xo x (<l = ®,E) 
thus extends from E — to the size of the sub-band A, approximatively E = 1.5 is 




(58) 



K xx (q,E)^2K°(q,E) 




(59) 



18 



0.8 



a- 0.4 









FIG. 2: DOS (CPA) for U=3.0, c=0.35 and \i = 0.02. 



this case. Transitions can also occur at higher energies between the band A and the 
impurity band B. A second peak is thus created that starts at an energy equivalent to the 
distance between band B and the chemical potential fi and ends at an energy equivalent 
to the distance between the higher edge of B and the lower edge of A. The result in 
this case is shown as the dotted line in Fig |3(c)[ U is large enough in this case for the 
large U limit (}28|) to be applicable. From the last two terms in (}28|) we can estimate the 
weight of the A-A transitions compared to the weight of the A-B transactions. Since the 
main sub-band is only half-filled, the ratio A-A/A-B is approximatively 0.9. Examples 
for smaller U are given in Fig |3(ajl and |3(b)1 (dotted lines). At U = 1.0 and U = 0.5 
(N e = 0.7), the two peaks can still be identified even though they are not separated by a gap. 



The effects of the vertex correction are depicted in FigO (solid line) for small, medium 
and large U, compared to the non corrected Tmxo x (q = 0,E) (dotted line). As predicted 
by ()29jl . the A-A transitions are slightly reduced while the A-B transitions are noticeably 
enhanced. This is to be understood from the perspective of the spin up and down 
propagating on the four-component alloy. As the spin up interacts with the down-spins and 
the spin down with the up-spins, it is to be expected that the n ± sites will play a major 
role in correlating the particles motions. For both the propagating spin up and down, the 
n sites are included in the impurity B sub-band, hence enhancing transitions of the A-B 
type. Similarly, the n° sites, with which no spin interacts, are included in the A sub-band 
of each species. Compared to the free particle cases, it is thus expected that the inclusions 
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(c) U = 3.0, n± = 0.0 and n = 0.02 

FIG. 3: Imaginary part of x xx (li E) and yp xx at q = and for N e = 0.7. Dotted line: Xq x without 
the inclusion of the vertex correction, solid line: Xo x with the inclusion of the vertex correction, 
circles (when appropriate): \ xx including both the vertex correction and the residual interaction. 

of correlations will reduce the transitions from one A sub-band to the other. 

The result of the approximation (JoUj) for \ XX {E) at q = that includes the residual 
interaction is plotted (circles) in its range of validity: Fig|3(a)l for U = 0.5, N e = 0.7 and 
n ± = 0.04, and for a split-band case at U — 3, N e — 0.7 and ri^ = 0.0, Fig |3(c)[ Whereas 
the inclusion of the vertex corrections generally reduces the A-A peak and enhances the 
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A-B peak, the correction ()59|) has different effect at small and large U. At small U the A-A 
transitions are enhanced compared to the value including only the vertex correction. At 
large U the A-A transitions have come back to their vertex correction value. In both cases 
though, the A-B transitions are strongly reduced. The approximation (|59|) is not valid at 
intermediate U. It in fact causes the imaginary part of x xx (E) to change sign between 
the two bands, hence producing unphysical results in this range. Therefore no results for 
X XX (E) were plotted in Fig |3(b)| 

An example of the singly and doubly weighted Green's functions obtained for the 
susceptibility is given in FigsE] and El for U = 0.5 and U = 2.0. In these pictures, 
K°_ + := K°, K 0+ := K 0n+ etc. We recall that x° xx — in the paramagnetic 

case. By construction, K° = K 0+ + K°~ + K°° + K 0± . Whereas x 0xx is an observable 
quantity and is therefore positive definite (and accordingly, the imaginary part of K° is 
always negative), nothing can be said about the weighted quantities. As can be seen in 
FigsE] and some weighted Green's functions contribute negatively for a part or for the 
whole spectrum. The n Green's functions contribute negatively to the A-B peak at all 
strengths of U. At low U, K 0+ is negative in the A-A band, and positive in the A-B 
band. At large U, the contribution from the n + quantities becomes comparatively small. 
For both the small U and large U case, the A-A peak is made up mostly of K°° and 
K 0o °, that is, from the n contributions. Similarly, the A-B peak is constituted mainly of 
the K°~ contributions. K 0± and K 0±± , due to the size of n , do not contribute significantly. 



Whereas the picture given by the singly weighted Green's functions for x~+ an d X+- 
corresponds well to the physical reality, the weighted residual interaction in the paramagnetic 
case can only be considered for Xxx- The limit m — >• was not obtained in a general case. 
Instabilities in the ground state, if they can occur, would happen when K dd (q = 0, E = 0) = 
—n + n~/U, with d G {+, — }. K dd (q = 0,E = 0) is going to zero when U increases. In the 
large U limit, we find 

This leads us to think that K dd (q = 0, E = 0) goes to zero with 1/U 2 , so that in agreement 
with the one-particle result, no instabilities should be found in the ground state. 
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FIG. 4: Imaginary part of the singly- and doubly- weighted Green's functions compared to the 
non-weighted case for U = 0.5, = = 0.35, n± = 0.04, \i = -0.05. 

The approximation (jo^j) is only valid in both extreme U limits and cannot be interpolated 
at intermediate U. However, it can be predicted that the peak at the energy of the maximum 
of the A-B transition will be reduced from its value obtained from the vertex correction. This 
is an expected behaviour considering the way we split the equation of motion. In the first part 
of the problem, when considering the propagation of two particles on a four- component alloy, 
the inclusion of the vertex correction led to the enhancement of the interband transitions. 
Due to the nature of the problem, it is indeed expected for the particles to become correlated 
when they meet on an site, the type of site that represents an impurity for both of them. 
However, the Pauli principle forbids the j and j spin to flip on such a site. When considering 
one-particle Green's functions, the Pauli principle is taken into account through the Fermi 
statistics. The vertex correction is obtained from the point of view of two particles moving 
in a lattice of potentials, rather than an up and down spins moving among other up and 
down spins. It is thus not clear whether the inclusion of the vertex correction conserves the 
Pauli principle or not. The inclusion of the last term in the equation of motion in the form 
of a weighted scattering correction appears to adjust the results of the vertex correction so 
that they comply to the Pauli principle. 
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FIG. 5: Imaginary part of the singly- and doubly- weighted Green's functions compared to the 
non-weighted case for U = 2.0, N 1 * = = 0.35, n± = 0.0, n = 0.02. 

IV. TWO-PARTICLE GREEN'S FUNCTIONS IN THE t-CPA 

The r-CPA treatment for the one-particle Green's functions is extended to the two- 
particle case, by using the two-particle CPA propagator up to a time r, thereafter the VCA 
propagator. The resulting Green's function is thus made of a CPA part that includes the 
vertex correction and of a VCA part, both weighted by factors depending on r. 

The two-particle Green's functions in the r-CPA can be built in a straight-forward manner 
from the one-particle case. The consistent way to build the two-particle Green's function 
(G^(t)) is to break it in a VCA part and CPA part, as for the one-particle case: 



For any time t, (G^(t)) = (G ] (t)G^(t)). In the one-particle case, the expression for the 
Green's function in time could be easily obtained in the VCA case. In the CPA 
reasonable form in time could only be assumed. As is discussed in I, this is due to the 
nature of the CPA poles. To obtain the two-particle function of the type (G^ 2 \E)) = 




(61) 
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du(G^(u, E-cu)), where G^(u, E - uj) = G\u)G l (E - u), one has to evaluate 

/oo 
dte iEt (G {2 \t)) (62) 
-oo 



oo 



DC 



rfte^(G( 2 ) CPA (t))+ / dte iEt {G^ VCA (t)) 



=:{G^ CPA {E)) =:{G(V VCA (E)) 

The VCA part of the two-particle Green's function is easily obtained. Using the definition 
from I 



J[E-e k -N- a U)T 



g«vcA {E)= t dte iEt Gf CA (t) = ^ — — (63) 

k K ' 2tt J t k K ' E - ejfe - N- a U 

the VCA part can be expressed as the convolution of the one-particle functions: 

/OO 1 
4ecw ca (£ (64) 
-oo JV 9 

and 0( 2 ) VCj4 (i) = £t( VCA )(t)£J- ycL4 (t). The factor for matching the two Green's functions at 
time r is included as usual in the VCA part, so that we have 

q 

with, from I, the one-particle VCA part that includes the normalisation f3k 

1 / IdE ct (£;)t \ 

<*™c*> == £ ( .-.-^ ) < 66 > 

The same approximation will be assumed to hold also for the CPA part. In that case, the 
procedure must be adapted to treat the vertex corrections correctly. The r-CPA says that 
the interaction is treated in the CPA up to a time r, from when it is treated in the VCA. 
In both cases, the interaction is taken in full. Therefore the vertex correction term must 
involve the full CPA Green's functions and not only the CPA equivalent of 



1 / 1 _ p ^(e)t \ 
«™==sfcb*E>) <67) 
The pure CPA two-particle Green's function is given by (|16|). A similar treatment to that 
for the one-particle Green's functions, where the energy dependence of E was ignored tem- 
porarily while the time integration was carried over, is applied here as well. The vertex 
correction to the free particle case, the denominator of ([160. is kept unchanged in the pure 
CPA form while the numerator becomes, as for the VCA case, a convolution on the r-CPA 
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one-particle Green's functions. The following approximation will thus be used for the CPA 
part of the two-particle Green's function in the r-CPA: 



(G k ^ CPA (e u e 2 )) 

Tf S G^ PA {ei)G q (e 2 



l-A( ei ,e 2 )^E 9 ^( ei )Gr A (e 2 )) 
A(ei, e 2 ) is calculated using the self-energy obtained from the full CPA equation. 



V. DYNAMICAL SUSCEPTIBILITY IN THEt-CPA 

In this section we put together the elements developed in the prior sections. The two- 
particle Green's function calculated so far does not take into account the extra [/-term of 
the equation of motion included in subsection IIIID1 In order to use the formalism developed 
above to calculate the susceptibility, the residual interaction between the two spins must be 
included in both the VCA and CPA two-particle parts. 



A. Construction 



The transverse susceptibility, without the residual term, is given by the two-particle 
Green's function K®(E), built from the disorder average of (|15|) . The approximation devel- 
oped so far gives 

K otCpa {e) = Kl VCA (E) + K° k CPA (E) (69) 

where K% (E) and K® A (E) are given by the convolutions (fTo^) of the two-particle 
Green's functions fj53j) and (p)%j) . 

The inclusion of the residual [/-term must be considered now. The correction for a pure 
VCA system is given by (f4"2*|) whereas the pure CPA requires the weighted result (|5Tj) . Once 
again, a straight-forward treatment would be to Fourier transform ()42|) and ()51|). then to 
integrate the CPA part from to r and the VCA part from r to oo. This would affect 
both the numerator and the denominator of ()42|) and (Jo"Tj) . The denominators contain the 
correction due to the interaction. It was argued in the context of the vertex correction for 
(G™ A (E)) that the correction is kept unchanged, since the interaction must be included 
at all time. An identical procedure is applied here, so that the time integration is carried 
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only on the numerators: 



K, 



tCPA 



(E) 



K 0VCA {E) 



+ 



(70) 



l + UKl VCA {E) 1 + K^ PA (E) 

K k VCA (E) is given by the convolutions (JT5j) of the complete two-particle Green's functions 
Gf VCA {e u e 2 ) = N-^ q G\lT (^i)Gf c ^(e 2 ). Kg PA (E!) and K*™^ are defined be- 
low. In order to simplify the notation, the index CPA will be abandoned for the weighted 



Green's function, and n + becomes just + (with identical definitions for the other species of 
the alloy). All functions are evaluated at the energy E. 



•"-life 



u 

m \ n 

U 2 
m 2 n + n~ 



rr 



(71) 



m \ n 



n + 



(72) 



K. 



OCPA 
2k 



K 



OCPA 



1 + K% PA 



U 
m 



n 



x 



+ 



U 2 



m 2 n + n 



+ [K k -YK2 + - 



2K° k + K k -K° k " 



(73) 



This formulation is not completely consistent with the rest of the picture. The denominators 
of both VCA and CPA parts do not yet contain the full interaction. We recall that (JT5j) is an 
integral from to the chemical potential /i. /i is obtained consistently from the one-particle 
r-CPA Green's function G rCPA (see I) and not from the pure CPA or the pure VCA. The 
problem can be retraced back to the development ()43jl relating K to K°. We note that 
the expression ()43|) is also valid in the case of the pure VCA, since in this case the factor 
(n-ni ~ n n^) is simply averaged to m, consistently with the VCA. Going back to the expansion 
(|43jl. it is clear that the term m in U /m^2 n (n^ — 7%f) comes from the inhomogeneous term 
2n5(t) (6jk{ctcii) — Sn(c^+Cji)) in the time version of the equation of motion (JHJ). The two 



26 



correlation functions involved are thus to be calculated at the time t = 0. If a pure VCA 

(CPA) system is considered, the m generated by the inhomogeneous term will be called 

m VCA ^ m CPAy 



m VCA _ N IVCA _ N WCA 



de 



de 



l)^J2MGr A (e)-G^ A (e)} (74) 



m CPA = N ICPA _ N ]CPA 



l)^T,MGi CP %)-Gr%)} (75) 



G%. VCA is given by the full VCA one-particle result and by the full CPA (as detailed in I). 
Hence, the expansion (}42|) for the VCA case should be written as 



K VCA (n p\ - K 0VCA (q,E) 

K iq > E) -l + U^K°vcA^ E) (76) 

In a pure VCA case m is simply equivalent to m VCA . The CPA case (15 1|) must also be 
re-casted in the same way, by subsisting m — > m CPA in (|71H7Hj) . If the case of a pure CPA is 
considered, m = m CPA . The approximation for the two-particle Green's function, involving 
the CPA vertex corrections and the extra [/-term of the equation of motion, is thus 

Kf PA (E) = 
K VCA {E) K^ A [m^ A ](E) 



l + U^K° VCA (E) 1 + K™ PA [m° p *\(E!) 

The chemical potential /x, being obtained from the whole G tCPA , will in general lead to an 
m cPA an( j m vcA q U j^ e diff ere nt from the values obtained in pure systems. 



It is crucial to keep the full interaction, as it guarantees the appearance of the spin waves 
at E = and k = 0. As ferromagnetism sets in, collective spin wave excitations are expected 
to appear. The magnetisation at k = can assume any direction in the averaged system and 
the energy cost E for exciting spin waves must be zero. This is also found in the r-CPA. The 
spin wave spectrum is generated almost entirely from the VCA part. As will be illustrated 
in the following paragraph, the CPA part at large U is negligible at low energy E. But even 
at lower U we recall that the CPA part does not create instabilities at low energy, so only 
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FIG. 6: DOS P ^ tCPA (E) and p lrCPA (E) using a simple cubic band with U = 6, iV T = 0.1, N l = 0.6, 
t = 1, w = 1. A\, B\, Aii -^3 denote the CPA sub-bands, V\ and V2 the VCA bands. The shaded 
area marks the occupied states. 

the VCA part of (J77|) is considered for the moment. In a ferromagnetic ground state, the 
following limit is obtained: 

VCA 

lim K k WCA (E) - — (78) 

Putting ()78)1 into (f7H|) . the result for finite and A; — > is 

limi^(i?) - (79) 

For small but finite k, we expect the excitation energy E to increase. As in the Stoner 
model, the spin waves will finally merge with the continuum. 

The transverse susceptibility Xk + {E) can now be obtained as x~k + (E) = ~Kk CPA (E)- 
In ferromagnetic ground states xt~i^) * s a ^ so required and can be obtained similarly to 
X k + (E). The transitions involved in Xk~ + (E) an d xV~~(E) can be understood by looking 
at FigEl a typical example at r = 1 of a ferromagnetic equilibrium situation in the r-CPA. 
X k + (E) calculates the transitions from occupied down-spin states to free up-spin states, 
whereas xt~(E) is concerned with transitions from occupied up-spin states to free down- 
spin states. FigEl shows the densities of states obtained from G]j CPA (E) and G k TCPA (E) 
calculated with the one-particle Green's functions obtained in I. In this example, U = 6.0 
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and N e = 0.7. The equilibrium was found to be for = 0.1, N^- = 0.6 and fi = 0.76. As is 
typically the case for large U, the chemical potential falls in the VCA band of the majority 
spin DOS. The sub-bands in the DOS have been labelled with Ai(A 2 ) and Bi(B 2 ) for the 
CPA sub-bands, and Vi(I^) for the VCA contribution. There is no transition between a CPA 
band and a VCA band. Therefore, the transitions of lowest energy happen between V% and 
V2 in x~k + {E)- At lower U, the CPA and VCA sub-bands combine in a complex fashion, as we 
have seen in I, and the CPA contributes at low energy too, though it never diverges at E — 0. 

The paramagnetic limit of (|77jl is again difficult to obtain, and will not be treated here. 
The CPA part presents the same difficulties in that encountered in the pure CPA. 

As for the VCA part, m/m VCA is varying in a complicated and self-consistent way with the 
species concentrations and the self-energies. It is however expected that this quantity will 
remain constant in the paramagnetic regime. It is therefore not unreasonable to assume that 
m/m VCA in the paramagnetic ground state retains the value it had before the transition. 
This value can be obtain numerically for all filling. For N e = 0.7, we find that m/m VCA ~ 
0.8. In the rest of this chapter we will concentrate on the application of Ij77|) on ferromagnetic 
ground states. 



B. The n problem 

Numerical calculations of the susceptibility requires the evaluation of n ± . This value 
is an equilibrium property of the system along with the chemical potential \i for a given 
concentration of spins up and spins down. We recall that the weighted Green's function 
formalism allows for the determination of n in the CPA, as presented in section 1111 C 11 
through (jHSJ)- (ESj) is written using the up-spins functions, but could equally be rewritten 
using the down-spins functions. In an equilibrium situation both formulations lead to the 
same solution, since the CPA leads solely to paramagnetic ground states for N e < 1. In the 
r-CPA the system can be polarised and (|33|) cannot be used, n cannot be deduced from 
the one-particle properties alone under general circumstances. In a general case n^ 1 could 
be determined self-consistently, or rather n~, since its correlation function corresponds to 
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the Green's function that is wanted for \ + 

n = N~ — n" 



dE ( — i ) ]m((c+c iT ;c+Ci i ))j 



= J dE (-Lj TmKuu(E) (80) 
Kiin(E) depends itself on N + ,N~ and n 

N a = f dE (—) TmF a (E) (81) 



For a fixed number of particles, (JSTJl is sufficient to determine N + ,N~ and /i, but (JHUl 1ST]) 
also form a closed set of equations with n^ 1 and fi to be established self-consistently. This 
simple formulation hides a very long numerical process. For this reason n has not been 
obtained that way for the susceptibility calculations in the polarised regime. For most cases 
of interest, such as when the system is close or beyond the transition from paramagnetic to 
ferromagnetic state, the value of the interaction will be fairly large. It is clear from FigQ 
that it is reasonable to put to zero under those circumstances. A has therefore the form 
given in f!25l) and (J2HJ), with ri^ = 0, iV + = n + and iV~ = n~ , in the numerical examples. 



C. Numerical results and discussion 

In this section the susceptibility is calculated using the result f)77jl for r-CPA at r = 1 in 
a ferromagnetic ground state, for a simple cubic band. The most prominent features that 
are discussed are the spin wave spectrum and the effect of the weighted development of the 
residual interaction on the A-B transitions. 



1. Spin waves 



When the ground state of the system has become ferromagnetic, we expect to observe 
the formation of collective excitations starting at q = and E — 0. Since the part of 
the spectrum at low energy is largely dominated by the VCA part, it is interesting to see 
first how the spin waves behave in a pure VCA case for a simple cubic band. FigsJZI and 
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FIG. 7: Spin wave part of the spectrum in the ferromagnetic regime of the pure VCA for U = 6.0, 
= 0.0, = 0.7. The imaginary part of x xx {l^) is plotted for q varying along the (1,1,1) 
direction 
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FIG. 8: Spin wave part of the spectrum in the ferromagnetic regime of the pure VCA for U = 1.6, 
= 0.05, = 0.65 . The imaginary part of x xx (Q:E) I s plotted for q varying along the (1,1,1) 
direction 



|S]have been calculated using the full VCA and its equilibrium condition at U — 1.6 and 
U = 6.0 for A^ e = 0.7. x xx {q^E) was evaluated for q along the (1,1,1) direction, from 
to the edge of the Brillouin zone at Q. At large U, Fig0 the down-spin band is full and 
separated by a gap from the empty up-spin band. The spin wave spectrum is characterised 
by a sharp peak at q — 0. The peaks observed at finite q broaden up as they enter the 
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continuum of excitations. For the parameters of FiglZl U = 6.0, N e = 0.7, we can see 
that for a q equal to 3/5 Q of the lattice vector, the spin- wave is already well within the 
continuum. When the ferromagnetism is not saturated, the up-spin band overlaps with the 
down-spin one. Figgis for U = 1.6, iV T = 0.05, N l = 0.65 and // = 0.28. As the band 
for the down-spin is more than half-filled, its Fermi surface crosses the zone boundary, but 
not on the (1,1, Indirection. The Fermi surface of the minority spin retains a spherically 
symmetric shape. The transitions allowed in calculating x~ + are given by the respective 
positions of the two Fermi surfaces as their centres are shifted by q. For a simple cubic 
band the result is strongly dependent on the direction of q. Transitions at zero energy 
can occur when the two surfaces intersect. FigJS] shows the spin- waves progressing along 
the energy axis. At two-thirds of Q, the spin-wave has completely disappeared into the 
continuum. At q = 13/15 Q, the spin- wave has come down in energy, and though still 
broad, starts to emerge from the continuum. The phenomenon is of the same nature as the 
antiferromagnetic instability present in the mean-field phase diagram for the simple cubic 
crystal structure. The system in the r-CPA in this range of parameters show a tendency 
to become antiferromagnetic. The separation of the spin-wave from the continuum is not 
observed, for example, along the (1,0,0) direction for the same set of parameters. 

The system is also in a ferromagnetic ground state at U — 1.6 and N e = 0.7 in the 
r-CPA with r = 1. The equilibrium is found at N r = 0.27 and N l = 0.43. At U = 6.0 
and N e = 0.7, the system is not saturated, and = 0.1, = 0.6. As expected, the 
ferromagnetism is weaker than in the pure VCA case, but some general features remain, as 
can be seen in FigsEJand^J At U — 1.6, the spin-waves are also seen to have decreased 
in energy at q — Q, but in that case the excitation is very broad and the peak is rather a 
resonance, FiglTUl The r-CPA causes the general broadening of the spectrum that is well 
visible in FigEI and ^2 by forcing the VCA Green's function to match the CPA Green's 
function at time r. 

2. The effect of the CPA part 

It has been discussed that the CPA part does not generally play a substantial part in 
the susceptibility at low energy. This is due to the choice of r that gives more weight to 
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FIG. 9: Spin-wave part of the spectrum in the ferromagnetic regime of the r-CPA with r = 1, for 
U = 6.0, N 1 * = 0.1, = 0.6, jj, = 0.76. The imaginary part of x xx (q,E) is plotted for q varying 
along the (1,1,1) direction. 
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FIG. 10: Spin-wave part of the spectrum in the ferromagnetic regime of the r-CPA with r = 1, for 
U = 1.6, A T = 0.27, N l = 0.43, fi = 0.34. The imaginary part of X cx (q, E) is plotted for q varying 
along the (1,1,1) direction. 

the VCA part, and the resulting density of states with the chemical potential falling within 
the VCA part. It has also been established that the weighted correction incurred by the 
residual term in U causes the reduction of the CPA A-B peak in the paramagnetic limit. 
This result is also valid in the ferromagnetic case, to the extent that this peak is completely 
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FIG. 11: Imaginary part of yp xx {q,E) in the r-CPA with r = 1, including the vertex correction 
but not the residual interaction. U = 1.6, iV T = 0.27, N l = 0.43. 




FIG. 12: Continuation of FigCHl at lar ger energies illustrating the elimination of the CPA A-B 
peak. Tmx xx includes the residual interaction, r = 1, U = 1.6,A^ = 0.27, = 0.43. 

eliminated by the weighting in r. FiglTTI show x\ X i before the residual interaction is taken 
into account. The CPA A-B transition band is visible at around E — 1.6. In FigEl where 
the residual interaction has been included, the effect is no longer visible. The broadening 
with q at around E = 1.5 is due to the broadening of the VCA part. Calculations for U = 6.0 
and Ag = 0.7 show for example that the CPA A-B transitions centred at around E = 6.0 
have become almost zero for any q. Compared to the paramagnetic case, correlations are 
also generally reduced by the decreasing number of scatters common to the two species. 



34 



The A-B peak is thus comparatively small even before the residual interaction is taken into 
account. 

VI. CONCLUSION 

In this paper we have extended our new treatment of the magnetism in the Hubbard 
model using the r-CPA developed in paper I to consider the two-particle propagators and 
the general magnetic susceptibility of the system. This required initially the study of a 
more complicated system within the CPA, effectively a four- component alloy since electrons 
of different spins see different scattering centres, to obtain the vertex corrections. The study 
was then extended to obtain the one- and two-particle weighted Green's functions within 
the CPA approximations. Such weighted functions allow, in principle, a self-consistent cal- 
culations of n* 1 (the concentration of doubly-occupied sites). This is only satisfactory in the 
paramagnetic regime, which is the only stable solution of the pure CPA. However in the 
r-CPA we found it necessary to use the CPA solutions in the ferromagnetic regime where 
inconsistent values of n 1 * 1 were found by the method. A fully self-consistent method based on 
the two-particle functions proved to be prohibitively complex and was not pursued, because 
for the magnetic cases of interest U is relatively large, in which case ri* 1 is always very small 
and the difference can be neglected. 

A full treatment of the magnetic susceptibility requires the residual interaction between 
pairs of electrons with opposite spins on the same site to be included. It is well known that 
this gives rise to spin wave excitations even within an RPA treatment of the Stoner (VCA) 
approximation. Within the r-CPA this requires the proper use of the weighted two-particle 
propagators, in part to ensure that the Goldstone mode at q = E = is preserved. Some 
further approximations are necessary to make the calculation tractable, which then provides 
a satisfactory treatment of the susceptibility. In the paramagnetic case, when treated in the 
appropriate limit, it shows an enhancement of transitions between sub-bands of different 
spins and a corresponding reduction of those between same spin bands. In the ferromag- 
netic case an appropriate spin wave spectrum is found. 

Although the treatment via the r-CPA reported here is restricted to a single band in the 
simple cubic lattice with less than half-filling, it displays a more general theory of the mag- 
netic state and excitations within the Hubbard model than any other previously reported. 
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Extensions to other forms of magnetic order such as antiferromagnetism which are expected 
for bigger filling and more complex band structures are straightforward in principle though 
no doubt complex in evaluation. It should be possible to refine the approximations made 
within the general framework of the physical idea behind the r-CPA. 

Acknowledgments 

AU acknowledges the support of the Swiss National Science Foundation, the Berrow 
scholarship trust of Lincoln College and the ORS award scheme. 

APPENDIX A: SINGLE WEIGHT FOR TWO-PARTICLE GREEN'S FUNC- 
TIONS 

We consider the evaluation of (G* n °CG^ on ). n~ is an "impurity" site for the up-spin, 
but a "host" site for the down-spin propagating, so that (G^ n ~°CG lon ~) = (G^ io CG loh ). 
This in turn can be decomposed as 

(G^ io CG loh ) = (G rio CG l ) - (G^ io CG loi ) (Al) 

This is particular to the four-component alloy for which the two propagating particles each 
see a different lattice of impurities and hosts. The second part of (|A1|) is already known, 
so that only (G^ lo CG^) needs to be determined. Using the definition (|3*TJl . the perturbation 
expansion 

G = P + PVG = P + GVP = P + PVP + PVGVP (A2) 
and the T-matrix expansion 

G = G + GTG (A3) 

we have 

Q io = (p-ig + p-^GTG - 1) (A4) 
A convenient set of equivalences are obtained from the Dyson equation G = P + PSG: 

p^G = GP- 1 = EG + 1 (A5) 

Putting (IA4IA5j) into (|A"lj) gives 

(G^ CG l ) = -[E^G^CG 1 + (S T G T + l) TG 1 } (A6) 
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where Velicky's original expression for the vertex T = (T^G^CG^-T^) appears in the last 
term of the second line of (IA6I). The result was found to be 2 * 



r n ( Zl , z 2 ) = \n) A( Zl , z 2 ) (n\ G (2) (z u z 2 ) \n) (n\ (A7) 

where T = ^2 n T n , which can thus be substituted in (jA6j) . The diagonal element of (|A6|) is 
therefore 

(m\ (G^ io CG l ) \m) = ^[S T (m\ G^CG 1 \m) 

n) (n| G^ \m) (n\ (G^CG^) |n)] (A8) 

n 

The whole weighted Green's function can now be obtained by Fourier transforms. Defining 

Qf := e ~ ikRm H (G^CG 1 ) \m) (A9) 
in 

Gf )loo ° := e~ lkRm (m\ (G^ io CG l ) \m) (A10) 

in 

a (2)n-oon- J^e^Rm ^| ^ n ~°CG lm ~) \m) (All) 
m 

The value of G^' is already known, (fTBj) . The result for Q^ 1000 j s 

£ T a fc + A (EU fc + F l ) Of] (A12) 



The Fourier transforms and Ak were defined in (II 7|) and (|18|). Gathering the results (|A1|) . 
(j2U) and (lAT2|) . we have finally 

ajo.--- = of (*±*t - (A13) 

Similar calculations can be carried out to obtain (G^ n °CG Lon } from (G^ ho CG loi ), and 
(G^ 00 ^^" ) from (G T/lo CG io?t ). We find 

y(2)nW _ + AF T + AD 



Gr °° n = <%> [^-jf- - U2 ) < A 1 1 ) 

( ^(2)n e 'oon _ ^(2) 

H 1 — ^ ^+^7^ 

It is more difficult to calculate (G^ on CG^ n °), with the weights swapped inside, as 
T l = (T^G^P^CG^ 1 ) appears in the evaluation of (G^CG 1 ) instead of T. T l will 
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be evaluated in the context of the double weights in the next section. With the knowl- 
edge of T it can be shown that for C diagonal, (G ]oi CG l ) = (G Vto CG l ). Elliott and 
Schwabe found that in this case, (G^ m CG^ ) is also equivalent to (G^CG^-™), so that 
(G^ on CG^ n °) = (G^ n "CG^ 011 ). This result is general for the other single weighting too. 
Using the expression for the vertex correction (|25|). the weights in (J34)) and (|A13HA15|) can 
be rearranged and give the results detailed in (j3Tlj) . 

APPENDIX B: DOUBLE WEIGHT FOR TWO-PARTICLE GREEN'S FUNC- 
TIONS 

The case {G Vn ~ n ~ CG [n ~ n ~) is equivalent to (G^CG lhh ). The Green's function sought 
after can be expressed as 

(G^CG lhh ) = (G^CG 1 ) - (G^CG loi ) 

- (G Vtl CG lio ) + {G Vli CG lii ) (Bl) 

The last weighted Green's function on the right-hand side is already known, but the three 
others must be separately determined. Using the perturbation expansion ()A2)1 and the 
T- matrix expansion (|A3J) . 

G Tii = — (P^GP- 1 + P^GTGP- 1 - P- 1 - V) (B2) 

The development (jB2|) can be substituted into the second term in the right-hand side of 
(IB1|) . After having eliminating the contribution containing (T), and applying (|A5|) . (|B2|) 

becomes 

(G^CG lm ) = M 1 +K 1 (B3) 
Mi = --^(V^C^G 1 - ^-(V^CG^^G^- 1 (B4) 

K x = (E^P-^CG 1 + p-^&Tx&p- 1 ) (B5) 
U 6 

V x := (T^&P- 1 CG 1 T 1 ) (B6) 

where V u is the interaction matrix for the spin a. The terms involving the disorder matrix 
are gathered in M.\. As discussed by Elliott and Schwabe^, for single-site disorder the 
matrices V* and are diagonal, and all contributions from j) vanish for i ^ j 
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when (|B3|) is projecting onto the real space. The total contribution from Aii is an order 
of magnitude smaller in N than /Q and can thus be neglected. Unlike the singly-weighted 
case (jA6j) . the vertex correction does not appear as such in the second term of JC\ (jB5|) . In 
order to obtain from (jB5|) a self-consistent equation, we develop the average as was done 
by Velicky for the vertex correction, using the effective wave approximation results and the 
definitions of Q and Q (Eqs.(33) to (36) in Velicky^) and get 

r\ = YsiQl&p-'CG^Ql) := ^r ln (B?) 

nm n 

Following through the procedure detailed in RefO, we find 

r ln = (T n T G T p- WTi) + (TjG^x^Ti) 

- (T^T ln G^) (B8) 

From (|B5|) . G^T\ n G^ can be extracted as 

G r T ln G l = U 3 PldP - ^Z l G^CG l P (B9) 

which leads to 

T ln = (T n T G T p- WTi) + U^TlPKxPT^) 

- (T^T ln G l T^ - ^Z l (T^CG l PTt) (BIO) 

Under this form, introducing the site projection T m = \m) t m (m\ allows us to extract A from 

Tin 

T ln = |n) A (n| {(n| G^P^CG 1 \n) 
+U 3 (n\ P/CiP \n) - E T E j (n\ &CG l P \n)} (Bll) 

By replacing ()B11|) into the original equation for /Q, (|B5J) . it is possible to solve the equation 
for (P/CiP). The following Fourier transforms are defined: 

)C lk = u 3 J2 e ~ ikRm (m\ p}c iP\ m ) ( B12 ) 

rn 

a k = J2 e ~ lkRm H G ] CG l P \m) (B13) 

Along with the definitions (|17I18|) of a k and A k , the solution reads 

~ tf&a k + AA k (tfa k + F l _ tf&a k ) 

/Cl fc = T^AA k (B14) 
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The complete solution for JCi(k) = Yl m e %kRm ( m \ ^1 l m ) can finally be found, after some 
algebraic manipulation 

r „, „ 2 /£ T E' +AD\ /E' Ft \ 

K 'W =G K tr- )\u + m (Bl5) 

As Elliott and Schwabe commented 3 -, it is consistent with the assumption that Ad% is neg- 
ligible to ignore the k- independent term. In the above (jB15|) . if the k-independent term 
(A£) + S^S^) is neglected, it reduces to the more convenient form 

KlW=flS (5^)(^) (Bl.) 

Since A*i is neglected, rtTTTTSfl is also the value for (JHHj), G 2iio \ 

A very similar calculation yields the same result for G 2mo , G 2um = G 2mo . G 2uo ° can be 
obtained along the same lines too, giving 

Gr°° = G 2 (^±^J 2 (B17) 
Combining (jHZJ), (jEEl and (fRT7j) and into (JHU) gives 

•-•-»-■ - = ( ST+ c/ Afi - £ ' Ei c/ : Ag ) 2 (B18) 

So we conclude that the weight of the doubly-weighted two-particle Green's functions is the 
multiplication of the simple weights, provided that M. Y can be ignored and (jB16|) is correct. 
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